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Abstract 

The interaction between a shock wave (attached to a wedge) and small amplitude, three- 
dimensional disturbances of a uniform, supersonic, freestream flow are investigated. The 
paper extends the two-dimensional study of Duck et al 1 , through the use of vector potentials, 
which render the problem tractable by the same techniques as in the two-dimensional case, in 
particular by expansion of the solution by means of a Fourier-Bessel series, in appropriately 
chosen coordinates. 

Results are presented for specific classes of freestream disturbances, and the study shows 
conclusively that the shock is stable to all classes of disturbances (i.e. time periodic per- 
turbations to the shock do not grow downstream), provided the flow downstream of the 
shock is supersonic (loosely corresponding to the weak shock solution). This is shown from 
our numerical results and also by asymptotic analysis of the Fourier-Bessel series, valid far 
downstream of the shock. 


’This research was partially supported by the National Aeronautics and Space Administration under 
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1 Introduction 


The interaction of freestream disturbances with shock waves is an important practical prob- 
lem, with particular regard to the operation of high-speed flight vehicles. Of particular 
interest are the consequences of these disturbances on boundary-layer receptivity, and this 
paper may be regarded as a step in this process. 

In the case of an isolated shock, Moore 2 , Ribner 9 and McKenzie & Westphal 4 showed 
that a single wavelength entropy, acoustic or vorticity wave, upstream of the shock, produces 
a disturbance comprising a mixture of all three classes downstream of the shock. When the 
shock is associated with the flow past a rigid body, the situation is much more complicated, 
due to the reflection and refraction of disturbances between the body surface and the shock 
itself, the latter distorting as a result, a process which further complicates the flow structure. 
This process was considered by Carrier 5 and Van Dyke 6 with particular regard to the 
problem of supersonic flow past a wedge performing small amplitude oscillations, the shock 
remaining attached to the wedge tip. These problems raise questions regarding the stability 
of the shock, and this aspect has been considered in the two-dimensional context (associated 
with wedge flows, the shock remaining attached to the wedge tip at all times) by Levinson 7 , 
Carrier 8 , Henderson & Atkinson 9 , Rusanov & Sharakshannae 10 and Salas & Morgan 11 . The 
overall conclusion is that if the flow downstream of the shock is subsonic (loosely classified 
as the strong shock solution), then the shock is unstable, in so far as disturbances grow 
downstream. If, on the other hand, the flow behind the shock is supersonic (loosely classified 
as the weak shock solution), then the shock is stable, with disturbances not growing in 
amplitude downstream. 

This aspect was considered in some detail by Duck et al. \ again for the case of super- 
sonic flow past wedges, subject to two-dimensional disturbances. This latter paper (where 
there is also a more detailed summary of previous work in this area) showed that two- 
dimensional disturbances, impinging upon the weak shock solution, do not grow downstream, 
but rather produce constant amplitude or weakly decaying oscillatory waves downstream. 
The conditions which distinguish between the two types of behaviours were determined and 
a decomposition of the far-field, suggested by the analytic solution, was made. The result 
was that there is a shadow region dependent on the incoming disturbance. If the direc- 
tion of propagation of the plane-wave acoustic disturbance generated at the shock by the 
plane-wave incoming disturbance intersects the wedge surface, then there is a non-decaying 
oscillatory pressure disturbance on the wedge surface. In any case, there are two weakly 
decaying single-wavelength oscillatory disturbances on the wedge surface that are due to 
the requirement that the shock remain attached. One is traveling at the mean-flow speed 
plus the speed of sound, and the other travels at the mean-flow speed minus the speed of 
sound. Thus, these are clearly regular acoustic waves which emanate from a point source 
(the wedge apex) of fixed temporal frequency as prescribed by the incoming disturbance. In 
addition to these more obvious effects, there is also a local field effect which is broad based 
in wavelength. While the first three modes of pressure disturbances at the wedge surface 
can indeed interact with the boundary layer, it is the last disturbance with the inclusion of 
short wavelength disturbances that should have the most effect on the growing boundary 
layer. In this latter paper, the side issue of shock stability to two-dimensional disturbances 
was also addressed. The aim of the present paper is to extend the ideas of Duck et al. 1 to 
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three-dimensional disturbances; indeed, little three-dimensional work of this type appears to 
have been undertaken in the past. 


2 Formulation 

Throughout this paper, we use subscript 1 to denote quantities upstream of the shock, and 
subscript 2 to denote downstream quantities. The wedge is taken to make an angle 0 with 
respect to the oncoming flow, with (x*,y*,z*) coordinates parallel and perpendicular to the 
upstream flow, with z* being the “crossflow” direction. The upstream basic flow velocity 
has magnitude U*, with Mach number Mi and density p{, and temperature T{. The ratio 
of specific heats 7 is assumed to be constant, and throughout this paper all our numerical 
results take 7 = 1.4. The velocity vector is written as U*(u , v, w), with respect to ( x *, y*, z *) 
coordinates, and the density is then written p\p : pressure as p\R*T{p (where R* denotes the 
gas constant) and the temperature is 7\*T. 

The Rankine-Hugoniot relations provide a link between conditions upstream and down- 
stream of the shock. For the basic (steady) flow, when the downstream flow is uniform and 
parallel to the wedge surface, the following classical result is obtained (see, for example, 
Liepmann and Roshko 12 ): 

tan(/? - 0) _ U 2 _ Pi _ (7 - l)M 2 sin 2 /? + 2 
tan/? Hi p 2 (7 + l)M 2 sin 2 /? 

where /3 is the angle between the shock and the wedge centerline (i.e., y* = 0) and Hi and 
« 2 are the non-dimensional velocity components perpendicular to the shock. 

The relationship (2.1) yields two possible values for /?, for a given value of 0 (for 9 < 
^max(-Mi)) with the so-called “weak shock” solution generally being characterized by super- 
sonic flow downstream of the shock, whilst the other “strong shock” solution is generally 
characterized by downstream subsonic flow. Note, however that there does exist a small, 
weak shock regime, close to 0 = 0 max where the downstream flow is subsonic, i.e. the 
downstream sonic line does not quite coincide with 9 max- 

We shall be concerned with the effect of small amplitude disturbances, which are intro- 
duced into the flow ahead of the shock. The corresponding two-dimensional results, as stud- 
ied by McKenzie and Westphal 4 , for example, may be extended into the three-dimensional 
contact, in quite a straightforward manner, when the z* variation is taken to be periodic (al- 
though other, more general z* variations could be accomplished using Fourier Transforms). 
We take e (< 1) to be a measure of the amplitude of the freestream disturbance, and it 
then turns out that the upstream disturbance may be classified into three distinct classes, 
just as in the two-dimensional case. 

(i) Acoustic waves: these are characterized by having a pressure perturbation, with 
corresponding perturbations in velocity, density and temperature, but no change in entropy 
or vorticity, to 0(e). Disturbances of this class take the form 

p = 1 + E + 0(s 2 ), (2.2) 
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u = 1 J, l£E + 0(e 2 ), 

7Mj 2 (oi + u> ) 

(2.3) 

v= ,+0(e 2 ), 

7M 1 2 (ai+u;) 

(2.4) 

w = - + 0(g 2 ), 

7M 2 (oi -1- uj) 

(2.5) 

p— Id — ~ + 0(e 2 ), 
7 

(2.6) 

and 

T= 1 + 0(e 2 ), 

7 

(2.7) 

where E is the normal mode exponential, i.e. 


E = expfioiXi + ia 2 yi + ia^z + icot}, 

with 

U) = Ox ± ^[ a \ + a 2 + a 3] 1/2 

(2.8) 

(2.9) 


being the frequency for the given wavenumbers. In the above, X\ and y\ are parallel 
and perpendicular to the upstream flow respectively (both perpendicular to z ); the non- 
dimensionalisation is carried out using one of the wavelengths of the disturbances as the 
typical lengthscale, e.g. by setting Oi to unity. The modes with the positive sign above 
in (2.9) are usually referred to as the slow modes, whilst those with the negative sign are 
usually referred to as the fast modes. 

(ii) Vorticity waves: these are characterized by having no density, temperature or pressure 
disturbances to 0(e), and thus the upstream flow takes the form 


and 

with 


V = £ 


u = l+eQ z , 2 , E + 0(e 2 ), (2.10) 

of + of 

+ (2.11) 

of + of of + ofj 

w = eQ x 2 Q2 -^ E + 0(e 2 ), (2.12) 

of + of 

p,T,p=l + 0(e 2 ), (2.13) 

u — — Oi- (2-14) 


The exponential term E is unaltered from (2.8), and the parameters Q z and must be 
specified. In comparison with the two-dimensional case, the three-dimensional case allows 
for an extra vorticity mode. 


3 



(iii) Entropy waves: these have no disturbance pressure or velocity components to 0(e ) , 
and may therefore be written in the form 

P — 1 + zE + 0(e 2 ), (2-15) 

T = l-eE + 0(e 2 ), (2.16) 

u,p = l + 0(e 2 ), (2.17) 

and 

v,w = 0(e 2 ), (2.18) 

with u defined by (2.14) and E by (2.8), again. 

In the following section of the paper, we consider the effects that these waves have on a 
shock wave attached to a wedge. The complication introduced by the presence of the wedge 
is that a single mode of the above type, upstream of the shock, will trigger modes of all three 
types, of all wavenumbers behind the shock; the only quantities preserved across the shock 
(provided e < 1) are the frequency parameter oj and the spanwise wavenumber 0:3. 


3 The Downstream Solution 

The flow downstream of the shock comprises the uniform flow solution plus a small ampli- 
tude (0(e)) perturbation, triggered by the freestream disturbances described in the previous 
section. We take non-dimensional coordinates parallel and perpendicular to the wedge (per- 
pendicular to the z direction in both cases) as x 2 and y 2 , respectively, with corresponding 
velocity components (u 2 ,v 2 ). 

We now write the solution downstream of the shock in the form 


«2 = U 2 + eu + 0(e 2 ), ( 3 . 1 ) 

v 2 = ev + 0 (e 2 ), (3.2) 

w=ew + 0 (e 2 ), (3.3) 

p = p 2 +ep + 0(e 2 ), (3.4) 

P = P 2 + £P + 0(£ 2 )i (3.5) 

and 

T = T 2 + eT + 0(£ 2 ). (3.6) 


Substitution of these expansions into the governing (inviscid) equations of motion, continuity 
and energy equations and equations of state then leads to 


Pt + U 2 Px2 + P2^X2 + P2Vy 2 + p?w z — 0, 
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PX2 = 0 , 


(3.7) 


Pv{ut + U 2 u X2 } + 


(3.8) 



(3.9) 


p2{v t + U 2 V Xi } + ^^2 Py2 — 

P2{Wt + U 2 U) X2 } + ~^j2 P z = (3.10) 

p 2 {T t + U 2 f X2 } - 1 -^ L {Pt + U 2 p X2 } = o, (3.11) 

and 

p = p 2 T + T 2 p. (3.12) 

Much of the success of the approach used by Carrier 5 , Van Dyke 6 and indeed by Duck 
et al. 1 was due to the ability of being able to split the solution into two components, the one 
including a velocity potential (which represented the acoustic wave component of the flow), 
and the other a stream function (which represented the vorticity wave component of the 
flow). It is clearly not possible to use the latter in the context of three-dimensional flows, 
however an alternative is the introduction of a vector potential. This concept, originally 
due to Poincare 13 , has been used in a number of fluid mechanics investigations over the 
years, although the appropriate boundary conditions have been the subject of discussion, as 
detailed by Aziz and Heliums 14 . 

Specifically we write the perturbation velocity vector u = (u, v , w ) in the form 

u = V0 + V A E, (3.13) 

where <j) is the velocity potential, and E the vector potential. The above represents a non- 
unique representation for u, and as such an additional relationship may be (arbitrarily) 
specified. The most popular choice, and the one that we pursue here, is that the vector 
potential can be required to be solenoidal, i.e. 

V • E = 0. (3.14) 

T his leads to the vorticity vector then being merely the quantity — V 2 E. Substitution of 
(3.13) into (3.8) - (3.10) leads to the equations 

V 2 |* + W„ + ^j|^J=0, (3.15) 

and 

V 2 {E* + U 2 E X2 } = 0. (3.16) 

It then follows that p satisfies 

P = —It, + (3.17) 

a s 

and elimination of the dependent variables p and T gives the equation for <p 

V 2 <t> = — ^[4>tt + 2 U 2 4> x t + U%(j) x 2l2 ], (3.18) 

a 2 
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where 


( 3 . 19 ) 


If we write E = (E^\ E&\ E^), then on account of our comments above, we impose the 
condition 

£% + E£ } + E™ = 0 . ( 3 . 20 ) 

A further quantity must also be introduced, namely xp'( y 2 , 2, t), which describes the displace- 
ment of the shock wave from its undisturbed state. 

This type of formulation then allows us to follow closely the two-dimensional approach 
adopted by Duck et al. It is now possible to write the general solution for 0, ip', and the 
in the following form, assuming boundedness at the apex of the wedge: 


i w t-'-^21X^-ia 3Z 00 

<t> = e u 2 ~ a ’ ^{a„cosh(i/0 2 ) + b v sm\i(vd 2 )}J v {k e r), 


( 3 . 21 ) 


. icjxo iwAnoaf . . 00 

E (») = e U2 v i (U i *i> 3Z y 


c^Mkeiy 2 ), 


( 3 . 22 ) 


iut-'^2Z!.+ia 3 z « * - 

xp =e v 2 ]T d v J v {k£y 2 ), 


( 3 . 23 ) 


where we have written 


( 3 . 24 ) 


- U) 

02 uj - qg 

P O' 


( 3 . 25 ) 


( 3 . 26 ) 


k = * 2 + ^, 


( 3 . 27 ) 


A = cot(/? — 0), 

| = yfx 2 - /3 2 , 

r 2 = xl- p 2 yl 


( 3 . 28 ) 

( 3 . 29 ) 

( 3 . 30 ) 


tanh 0 2 = /?y 2 /x 2 . 

If we impose impermeability on the wedge surface, we immediately require 


( 3 . 31 ) 


K = 0 Vi/, 


( 3 . 32 ) 
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whilst if the shock remains attached at the tip, then d 0 = 0. As noted by Duck et al. 4 (for 
example), it is also reasonable to set the v = 0 coefficients of the other terms in the series 
to zero, which correspond to transient-type modes, and so we write 

flo = C™ = cf = cj 3) = 0. (3.33) 

The key results for the velocity components and pressure are then 


u = — e 


* * iU^u>Z 2 , „ OC 

tut- -f^f +ia 3 z 


£ 

v=0 


2ia>f7 2 

k e {3 2 a s 


a„ 4- o [/ _i(l + <5 v _x,o) — G 1/+1 


cosh v>6 2 Jv(k e r) 


+€ 




£ 

i/=0 


2i\uia s 

le& 


4 3) + i (<SL - 4_i(i + 5,-1, o)) ) - *M 2) 


(3.34) 


iut-iSgEf+ia,, ^ k e 0 


v = e i 1 


53 -7r(0v - 1 + o«/+t) sinh(i/0 2 ) J„(fc e r) 

i/=0 ^ 


5: ^ 


iu/xo iutXya. 


T±T,+ ia * z 


IU 


z=0 


Ut 


Ji/{ke^y2) ? 


(3.35) 


. , _ OO 

iwt— g 


it) = e ^ a * 


i^duJv{K r ) c °sh ^#2 


iu/ip tuJAy2 a ^ OO 


— e 


and 


iurt+ias*--^- 


>E 

t/=0 


t/=0 

j2) j ^ 2 zAu^Q (5 

^ + 2 \ fc e E/ 2 /3 2 


4 1J + | (c^x — c[ / 1 3 1 (l + <5^_i,o)) 


(i 


Ji/(.k e ^y2 ), 

.36) 


p=— — e ^ s 2 - 

za s [,=0 


2wt> U 2 , : , U 2 

- T- a, , H + </i/-i,o) a ,+i 

k e 3 2 a s a s 


cosh(v'Q 2 )J„{ker), 

(3.37) 

where a_i = c ( ”i = eLi = 0, and where (5 n>m is the Kronecker delta. 

Although the solutions above satisfy the impermeability and apex conditions, as yet 
these solutions do not satisfy the unsteady Rankine-Hugoniot conditions (McKenzie and 
Westphal 4 ) on the shock; following (and extending) the results of Duck et al. l , these take 
the form . . 

r / 1 \ 1 

(3.38) 

(3.39) 

(3.40) 


u' 2 + tt 2 — = 1 R 1 4 {UiPl+u\), 


P 2 


P 2, 


P2 


, 1 _ p2 , a 2 s P2 Pi 

u 2 + n u 2 H 7~z = + -UlPl + „ , , 2 _ 5 

2 2 P2 2yu 2 p 2 2 2'yMfui 


u 


a 2 p2 al P2 


2 (7 — l)u 2 p 2 ' (7 -l) a 2 P 2 

* + 


V u 2 J u 2 


P 1 


(7 — l)M 2 u 2 (7 — l)M 1 2 u 2 ’ 
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and 


(3.41) 


with 


<4 = + (Ui - tZ 2 )l/£, 

«4 = + («i - u 2 )V>' , (3.42) 

= ip' t + Viifc. (3.43) 


In equations (3.38)-(3.43), a denotes the unperturbed distance along the shock, (u',v',w') 
denote the velocity perturbations perpendicular to the undisturbed shock, parallel to the 
shock and perpendicular to the z direction, and parallel to the z direction, respectively, and 
(u,u,0) denote the corresponding base flow velocity components. Eliminating p 2 from the 
above equations, and then utilizing (3.34)-(3.37), yields the following four equations for each 
of the v > 0: 


cos(0-0) j-y 


2iuiU 2 

Ue/3 2 a, 


a„ + a v - i(l + 5 v _j > o) — a„ + x 


cosh u9q 


+ 


l 2 ( ~ + { " ' a,(1 + *— >) - 


+ sin (0 - 9) + (Wi)sUi(4) + 0c« + 

— (Hi — Hi) sin(/l — 9)~ ( - - d u + £ (d„ + i — d v - i(l + 5 v _i,o))| = Rl l \ 


k e (5 2 a s 


UH- 


sin(/? - 9) 
k e ( 2i\iua 


2iujUo 


a, v + a„_i(l + ^-^o) — <Zi/+i 


2 V k e U 2 (3 2 
cos (/? — 9) 1 


2 [k e /3 2 a 

C? + i (cgjx - cSx(l + (Wo))) - ^cl 2) 


cosh u9 0 


(a *- 1 + (Wx) sinh(i/0 o ) + + 77- c£ 


.(3) 


{ 


+ > < vx sin(/? 9) — 


2iu>U 2 \ 

k e (3 2 a s 


U 2 


d n u + £ (d|/ + x — d 1/ _x(l + — x,o) ) 


(3.44) 


4- iu)d v 


+ 


\ A z \ \ k e \ 2 iu U 2 U 2 ) ) 

1 B 3 ) \Jj2 a - + — °-i( 1 +^-i,o) - -a„ + ,|cosh(^ 0 )| 

= I *!?’ \ 

\ Kl 3> ) ’ 


(3.45) 


and 


ia 3 cosh(^ 0 )a 1 ,-^c[ 2) -^(-|^^c' / 1) +| (c^x - cSx)(l + ^-i,o))-ia 3 (wi-u 2 R = R {4) . 


U 2 " 2 k e U 2 0 2 


(3.46) 
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These equations must also be augmented with the condition arising from the solenoidal 
condition on the vector potential namely 


I^(D 4. bl(. 

* o V 


2i\Coa 


— — c 


l h 


k e U 2 (3 


|c<?> + 4 (<£>, - <S)(1 + + «,<£» = 0. (3.47) 


Terms used in the above are 

tanh^o = j3tan(@ — 9), 


B 2 = 




B x = 1 + 


at 


(7-1)®!’ 


a* _ Hi" 
(7 - 1)«2 “2. 

D _ 7 fl S 

3 (7 - l)Ba 


(■ 



(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 

(3.54) 


The as in Duck et al. l , are to be determined from the freestream conditions. In 
particular in order to write the exponential terms, arising from the upstream solution, in 
terms of Bessel functions the following is particularly useful: 


exp[i(ai cot P + a 2 )(cos0 4- A sin 0 ) 2 / 2 ] = e u *~ a * 92 e *sm (3.55) 


where 



E 


t /=0 




4 - (-lye-** 

1 + fiv,0 


Jv (^e^!/2)> 


(9 = sin' 1 


(ai cot P + a 2 )(cos 6 + Asin0) + 

ki 


1 


(3.56) 


(3.57) 


Other details are routine, and omitted for reasons of brevity. 

The combined system (3.44)-(3.47) then represents a closed (recursive-type) system de- 
termining the a v+ i, cJSl d u+ 1 . 
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4 Some Numerical Results 


The first set of data we present is for the particular case Mi = 5, 9 = 25° (weak shock 
solution), for the particular case of a slow acoustic mode with «i = 1, a 2 = 0, and a 3 = 1. 
Fig. la shows results for the perturbation pressure on the wall ( y 2 = 0). Here, and in all 
cases the solidus denotes the real part of a function, the dashed line the imaginary part. 
This figure presents a picture reminiscent of results found in a number of corresponding two- 
dimensional cases by Duck et al. 1 , namely that of an oscillatory-type nature downstream. 
Fig. lb shows the corresponding perturbation shock location; this too takes on a similarly 
oscillatory nature. Figs. 2a, 2b show the corresponding results for the fast mode case (other 
data remaining the same as for figs. 1). These figures are to be compared with fig. 2 and 
fig. 3 of Duck et al. 1 where the same conditions apply except that a 3 = 0 in those graphs. 
The difference between the two-dimensional results and the three-dimensional results are in 
whether or not the disturbance produces an acoustic wave behind the shock that intersects 
with the surface. In response to this particular slow mode, the pressure response at the 
surface for the two-dimensional disturbance decays algebraically while the shock position 
shows a characteristic single wavelength response in the far-held limit. On the other hand, 
the surface pressure response to the three-dimensional slow-mode disturbance does not decay 
and shows a multi- wavelength pattern as is to be expected if a decomposition of the far held 
behaviour applies to the three-dimensional case as it does for the two-dimensional case. 
The shock position shows a two-wavelength structure as would be expected by the primary 
acoustic disturbance reflecting off the wedge surface and interacting with the shock from 
behind. There is very little difference between the responses to a two-dimensional fast- 
mode disturbance and a three-dimensional fast-mode disturbance. A search in parameter 
space would lead to the same conclusion as in the two-dimensional case: either the surface 
pressure has a non-decaying component and an algebraically decaying component or it has 
only an algebraically decaying component. This is signihcant in that the solution in the 
absence of a wedge indicates that the choices are between only a non-decaying response and 
an exponentially decaying response (see Hussaini et al. 15 , Jackson et al. 16 , Lasseigne & 
Hussaini 17 ). This issue was pursued in detail in Duck et al. 1 where it is determined that 
the condition that the shock remain attached to the apex of the wedge is the source for 
the algebraically decaying portion of the response - the apex is in effect a point source for 
additional acoustic waves. The interaction of these algebraically decaying acoustic waves 
with the shock in turn produces an algebraically decaying portion in both the entropy and 
vorticity modes. Another reason that it is important to determine that the pressure has 
only a decaying or sustained oscillatory behaviour downstream is that, as mentioned in the 
analysis of the response to the two-dimensional disturbances, the potential functions have 
components of exponential growth in them. The calculation of the physical quantities of 
pressure, velocity and vorticity are therefore dependent on some fortuitous cancellations, or 
they themselves would be exponentially growing. In the following section we investigate the 
downstream behaviour of the three-dimensional perturbation solution in some detail. 
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5 The Far Downstream Behaviour/Stability Problem 


The results of the previous section show that a general growth in physical quantities, as x 2 
increases does not occur even though the potential functions have an exponentially growing 
component in them. This is similar to the two-dimensional results as considered by Duck 
et al. l , in which it was shown by examining the recurrence relations that only decaying or 
sustained physical solution oscillations were possible. The same type of analysis is possible 
for the three-dimensional case, but indeed has its own peculiarities as will be seen. 

In order to analyze the x 2 — »• oo behaviour, we consider the homogeneous solution of the 
system (3.44)-(3.47), and consider the limit as v — > oo. In particular, we write 

4 n) « C5 n) /r, d v « DqK u , and a„ « A Q K v e~ u6 ° , (5.1) 


where K is an eigenvalue whose precise value (including location in complex space) serves 
to determine the ultimate downstream behaviour of the disturbance. If \K\ > 1, then using 
the generating function for Bessel functions, we have (for example) 


ES,i dMz) a Z? = - oo D 0 K‘'Uz) 

« D 0 e^ K ~^ z as \z\ oo. 


(5.2) 


Thus exponential growth occurs as \z\ — » oo if 9? (K — ^)} >0. If we substitute (5.1) 
into (3.44) - (3.47), allow u — » oo, and discard the inhomogeneous terms, then we obtain the 
following five linear homogeneous equations for the coefficients Aq , and Do'- 


\i0A> - ~Cf> - h 

A U2 


(- 


iXuJds If 1 






k e U 2 0 2 2 
cos ((3 — 0)\ — 

+ k e 


Cq 1) - ij3(ui - u 2 )D 0 = 0 , 

Cj 2) +*0Cf = 0, 


1 

\ 2l ^ + eh e -°. K ] 

l 4 

k e (3 2 a s K 


, iAuJCls 1; /Tr 1 . 


Ao 

C S 3) - i/JCf 1 ] 


+ sin(,3 - 9) + e~ 9o K)Ao + i0Cj, l> + ^C{ 3 > J 

-(6. - «j) sin(0 - 9)k . + §(* - ^)} A, = 0. 

( 


A\ 

E h 


sin {/3 - 0) 


U2MJ l + e*_ e ^ K \ Ao 

4 \k e 0 2 a s K ) 


(5.3) 

(5.4) 


(5.5) 
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~\~kf 


iXurn, 1> 1 

L k e u 2 p 2 ^ [K k\ 


dP-ifidP 


- cos(/? — 0) 


+ 


^2 

B 2 

+ 


^3 

B 3 


M<£ + e+K)* + iW + %0P 

( iuU 2 \ i 1 \ 

Ufc + 2 (K ~K ) ) JrlU ’ 

^_^VU) =0 . 

& S a s J J 


Vi sin (j3 - 9 0 ) 
f k e ( 2 id) 


k. 


D r 


2 \k e 0 2 


(5.6) 


This equation is nonlinear in K but can be converted to a tenth-order linear, generalized 
eigenvalue problem, which may be solved using the QZ algorithm by using the additional 
variables 

io = KA 0 , Ci n) = KC { f\ 

(5.7) 

D 0 = KDq. 


It was shown by Duck et al. 4 that in the two-dimensional case, for situations in which the 
downstream flow was supersonic, there were four imaginary eigenvalues (all with \K\ < 1), 
and two complex eigenvalues given by 


A- = 0ia»±[t/|-aJ]‘ /2 }. (5.8) 

However, although these complex eigenvalues are such that \K\ > 1 , in both cases D 0 = 0, 
and due to some “fortuitous” cancellations, all physical quantities remain bounded as x 2 
oo even though it is determined that the potential functions themselves grow exponentially. 

In the three-dimensional case, we may expect, ten eigenvalues instead of the six in the 
two-dimensional case. However it turns out generally that there are just eight distinct 
eigenvalues (two of the eigenvalues being double eigenvalues). Generally we find: (i) four 
imaginary eigenvalues, (ii) one complex conjugate pair of double eigenvalues of the form 
±K t + iKi , and (iii) one complex conjugate pair of eigenvalues also of the form ±K r + iK t . 

For 9 = 25°, Mi = 5 e*i = 1 and a 2 = 0, the four imaginary eigenvalues have magnitude 
less than unity even as a 3 increases from zero. The eigensolutions associated with these 
eigenvalues do not contribute to any growth in any of the quantities calculated. Most in- 
teresting, however, is the variation of these four eigenvalues as the wedge angle is increased. 
All four eigenvalues approach K = i at the maximum wedge angle. 

The (imaginary part of the) eigenvalues described by (i) are shown in fig. 3a, whilst the 
real and imaginary parts of the eigenvalues described by (ii) and (iii) are shown in figs 3b, 
3c respectively. Here we have just shown the eigenvalues with the positive value of K r (the 
other eigenvalue, corresponding to — K r + iK x may obviously be simply deduced). We next 
consider the effect of varying 9 , for the parameters Mi = 5, au = 1, a 2 = 0, and o 3 = 1. The 
variation of Ki of family (i) is shown in fig. 4a, and the variation of the real and imaginary 
parts of eigenvalues described by (ii) and (iii) is shown in figs. 4b and 4c respectively. Again, 
just the positive values of K r are shown. It is seen that the double eigenvalues described by 
(ii) are the only ones with \K\ > 1 and therefore lead to potential exponential growth. It 
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turns out that an analytic description of the eigenvalues described by (ii) is possible. The 
result is 


ika s 
k e 

It is also possible to obtain analytic results 
these eigenvalues, namely 



m - a ’ k2 


1 1 / 2 ' 


kl 


(5.9) 


for one set of eigen-coefficients corresponding to 


Q 


(i) 


{_°2kA( Ke -eo + ^) + a 3 ~ j<) ~ ^f} 


^ + g 


-U, 


iXka, 


U 2 


-¥(*- Tli' 


(5.10) 


cP = 


iVi 

UJ 


-ia^Ao + C 0 


\ka s k e if 1 W\ 

vr~~\ k ) Jr 
c «3 ) = ^{M(^4)^ +!Q3C -}. 


D o = 0. 


(5.11) 

( 5 . 12 ) 

(5.13) 


The above solution readily reduces to the two-dimensional solution in the limit (*3 — > 0. Some 
understanding of the origin of this other family of coefficients for this eigenvalue, and also 
of the additional complex eigenvalue (iii) may be made by considering the two-dimensional 
limit, (*3 — > 0 and seeking a solution to the system with Aq = Dq = Cq = 0. After some 
algebra, we find the following (four) eigenvalues: 


A ika s 
k£U 2 


±{ 


A ika s 

7Wh 



(5.14) 


together with 


^(1) _ U 2 kety f j/ 1 2A ika s 1 

Co {K ~k~ 


(5.15) 


Two of the above values for K correspond to (5.9), and this partly explains the origin of the 
double eigenvalues. As a 3 increases from zero, it is expected that A 0 , D 0 , and Cq will no 
longer be zero for this branch. The other pair of eigenvalues correspond to the eigenvalues 
described by (iii). 

Again it is important to note is that with the exception of the eigenvalues (ii), all eigenval- 
ues of K correspond to \K\ < 1, and thus are of limited physical significance. Interestingly at 
the maximum value of 0, (corresponding to the maximum wedge angle for attached flow at a 
particular freestream Mach number), all eigenvalues take on the same value, K = i. Turning 
now to the effects of the eigenvalues described by (ii), again, as in the two-dimensional case, 
it turns out that although seemingly these eigensolutions correspond to unstable (grow- 
ing) downstream solutions, in fact the evaluation of physical quantities leads to solutions 
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downstream which remain bounded. This is best illustrated by defining the following set of 
functions: 


*(») 


iut 


4> = A 0 e 


t wr 2 tu>Ay 2 <*, , • OO 

1/2 E ^J v {kiy 2 ), 

u——oo 

(5.16) 

iwC/ 2*2 , ... OO 

^ £; {Ke’^YUkr), 

(5.17) 


where K is defined by (5.9) and we write E = (E^\ E^, E^). It would be expected that 
E and <p determine the far-downstream behaviour of the flow which grows downstream. 
However inspection of the solution above, in the cases considered, leads to the conclusion 
that 

V0=-VAE. (5.18) 

This then leads to all velocity components (and hence all physical quantities) remaining 
bounded downstream. Numerically, it is determined that all physical quantities remained 
bounded no matter what the parameters selected for the calculation, thus the above must 
hold for all eigensolutions with \K\ > 1. 

In summary, therefore, we see that provided the flow downstream is supersonic then 
disturbances will not grow downstream. This boundary of downstream behaviours is inde- 
pendent of the crossflow wavenumber a 3 , and occurs along the ’line’ where (3 = 0 (where 
k and k e simultaneously become singular, and also change from being real to imaginary 
quantities). The result is that all of the Bessel functions in the series (3.21)-(3.23) will grow 
downstream if the downstream flow is subsonic (although we note the comments of Salas & 
Morgan 11 stating that strong shock solutions could perhaps be stable under a different set 
of boundary conditions) . 

We now address the nature of the downstream response to sustained excitation (i.e. 
acoustic waves, entropy waves, or vorticity waves). Arguments similar to those employed by 
Duck et al. 1 are again appropriate. The frequency of the disturbance downstream of the 
shock is determined by the frequency of the disturbance upstream of the shock and both 
wavenumbers of the disturbance tangential to the shock are also fixed by the deflection of the 
shock produced by the upstream disturbance. There is a finite range of tangential wavenum- 
bers that exclude the existence of plane acoustic waves which propagate downstream. If 
the tangential wavenumbers lie outside of this range, then a plane wave propagates at an 
angle $ p to the shock. In the case of the tangential wavenumbers being within the excluded 
range, an acoustic field is generated that decays algebraically rather than exponentially as 
would be the case in the absence of the wedge. If the tangential wavenumber is outside of 
the excluded range, then the solution in the presence of the wedge exhibits two types of be- 
haviour, dependent upon the angle of propagation of an acoustic disturbance in the absence 
of the wedge. If the angle of propagation Q p is greater than the angle between the shock and 
the wedge 13-9, then the pressure disturbances generated at the shock intersect with and 
reflect from the wedge surface, leading to a non-decaying pressure field everywhere between 
the shock and the wedge surface, with a superimposed algebraically decaying pressure field 
owing to the requirement of an attached shock and zero velocity normal to the wedge sur- 
face. If the angle of propagation 6 P is less than the angle between the shock and the wedge 
surface, then the pressure field is divided into two regions by a ray emanating from the apex 
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and parallel to the direction of propagation of the pressure disturbance in the absence of 
the wedge. Between the shock and this ray there is a non-decaying pressure field, with the 
aforementioned superimposed algebraically decaying component. Between the ray and the 
wedge surface the pressure field is algebraically decaying and there is no component with 
sustained oscillations. This was referred to as the shadow region by Ducket al. l . 

The various downstream limits are determined by the value of 9. If 6 is real, then the 
wavenumber of the imposed disturbances lies within the excluded range, and the pressure 
along any ray 9 2 = constant decays algebraically For 9 2 = 9 0 the solutions does not decay 
leading to the appropriate oscillatory behaviour necessary to match conditions at the shock. 
If 9 is complex, then the condition 


(oi cot 3 Q! 2 )(cos 9 + A sin 9 ) + 


> 1, 


(5.19) 


is satisfied. We also note that for equation (3.55) to be satisfied then sin 9 must remain real, 
which requires that 


9 — ±7r /2 — i9i , 


(5.20) 


where 


9i = cosh 1 { ± 


(qi cot 3 + a 2 ) ( cos # + A sin 9) + 

hi 


(5.21) 


91 is related to the angle between the direction of the shock and the direction of the acoustic 
disturbance 9 V through tanh 9 t = /3tan 9 p . If 6 t is greater than 0 o ,_then the generated or 
transmitted acoustic waves intersect the wedge surface. For 0 < 9i < 9$. there are two 
regions: a region of sustained oscillatory behaviour of the pressure field near the shock when 

9 2 > 9i , and a region of algebraic decay of the pressure field near the wedge surface when 
9 2 < 9i 

These features are common to those found in the two-dimensional case by Duck et al. 
F Additionally, we can expect that the downstream (forced) behaviour will take the same 
form as in the two dimensional case. Taking the wall pressure, for example, we may write 
the decomposition 


iujUox 




P = Poe 


\ 2 l j +tfc e cosh(*?,-0o)s2 


+ Pix 2 ^ 2 e 


( ,u, *2 


+ P 2 x 2 1//2 e 


f - iu>x .2 jZL) 


i.e., a solution with an oscillatory component, two decaying acoustic waves, and 
decaying component, Q(x 2 ). 


+ Q (% 2)1 
(5.22) 
a (faster) 


6 Conclusions 

The interaction between a shock wave and three-dimensional freestream disturbances has 
been considered, in particular the nature of the flow far downstream. The somewhat fortu- 
itous cancelations that were found to occur in the analogous two-dimensional work (Duck 
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et al. *), again intriguingly occur, yielding a non-growing physical solution downstream, 
thus confirming the stability of shocks with downstream supersonic flow, in line with the 
widely-held belief regarding the stability of such shock waves. Our results also point to the 
‘instability’ of shocks with downstream subsonic flow, since in that case the argument of 
the Bessel functions in (3.21)-(3.24) becomes imaginary, and as a result exponential growth 
will be expected downstream. Finally we note that the location of the boundary between 
growing and non-growing downstream waves is unaffected by the crossflow wavenumber a 3 . 
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Figure la. Wall pressure perturbations, M\ = 5, 6 = 25°, (weak-shock solution), slow 

acoustic mode, = l,e* 2 = 0,o 3 = 1. 
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Figure 2a. Wall pressure perturbations, = 5, 6 = 25°, (weak-shock solution), fast 

acoustic mode, c*i = 1 , q 2 = 0, 03 = 1. 
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